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Abstract 

The results of study of the gravitational fragmentation in the interstellar medium (ISM) by 
clump-clump collisions are presented. We suggest, that collision of clumps, that are subparts of 
Giant Molecular Clouds ( GMC) may be on of the basic mechanism, which result to ISM frag- 
mentation and define the dynamical as well as statistical characteristics (e.g. the mass spectra) 
of protostellar condensation. 

In the present paper, we describe our 3D SPH-modeling, in isothermal approximation, of 
supersonic collisions of two identical clumps with a few variants of initial impact parameters ((3), 
that cover the wide range. 

Our results shown, that at all (3 in system began intensive fragmentation. The resulting 
fragments mass function depend from initial impact parameter. The obtained mass spectra have 
the slopes in a good enough agreement with observational data for our Galaxy - especially for 
large impact parameters, which are more realistic as for large clumps ensembles. 

KEY WORDS ISM, fragmentation, molecular clouds collisions, clumps mass spectra, SPH 
method 



1 Introduction 

Nearly the half mass of gas in our Galaxy are concentrated in giants molecular clouds (GMC) with 
typical sizes in tens parsecs, masses 10 5 — 10 6 9Ji Q and temperatures T < 20K [H, [3, [3]. The mean 
number density in GMC is equal to 10 2 — 10 3 cm -3 , but in some regions it can reaches up to 10 6 cm -3 
[3]. There are about 4000 such objects in Galaxy and this objects mainly concentrate in the Galaxy 
disc U0, 3. 

The GMC tightly connect with some other objects in our Galaxy - such as IR-sources, masers, 
compact HII zones, T Tau stars, Herbig-Haro objects, O-association etc. This connection show the 
clear evidence of star formation process in GMC at present. 

The hight resolution radio observations of GMC in molecular lines, show their clear fragmentary 
structure on all hierarchy levels down to parsec's fractions. In first approximation we can describe 
the GMCs as a complex of smaller and denser condensations (clumps), embedded into tenuous gas 
medium. 
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The typical clump have a mass in 10 3 — 10 4 9JTq and size in a few parsecs. It content yet more 
small and dense cores with masses in several solar mass and sizes in parsec's fraction 0]. As a good 
examples of such clouds complexes we can brought the p Oph 0|, M17 SW [H] and Tauri molecular 
clouds complex 



The numerous observational data about molecular clouds indicate the complex gas flows and 
clumps chaotic motions inside GMC. There are strong turbulence inside clumps, but tenuous gas 
"background" have a more regular kinematics. The dependence of velocity dispersion a from size 
of consider region can be approximated by the power law a oc L 0,38 [3]. Here a is total 3D-rms 
velocity deviation for all internal motions including large-, small-scale and thermal gas motions 
(for H2 clumps at temperature in a few tens of Kelvins the mean quadratic velocity of thermal 
motions is negligible by comparison with two first components), a is about 1 km/s for R ~ 2 pc or 
Wl « 2000 9JT . 

This discovered dependence is very similar to Kolmogorov law for distribution of velocity dis- 
persion at subsonic turbulence (a oc L 0,33 ) and it is possibly, that all seen us interstellar movements 
is parts of a single hierarchy of interstellar turbulence. 

At that time the rms velocity dispersion of molecular clouds and clumps itself have mainly flat 
spectrum in wide mass range (up to four orders) with average value around 5 — 10 km/s @, 0]. 

Usually, the GMC, havin g m asses 10 5 — 10 6 9JIq, may contain a tens of clumps with typical 
masses by order 1O 3 -1O 4 9JT [lfl|. Of course, chaotically moving at supersonic velocities, they can 
collide each with other. Many authors consider this mechanism as one of the possible source of ISM 
fragmentation and formation of initial mass function (MF) of protostellar condensations. 

It is interesting to note, that as GMC, as clumps and cores inside its have almost the same mass 



spectra in mass range from 10 4 to 10 4 Off© 111 ]. This may point to common mechanism of their 



formation. For example, just such result was obtained in special N-body modeling of collisional 



buildup of molecular clouds [l2j, that gave observational mass- and velocity spectra for obtained 
objects. 

So, interclouds collisions have a great interest for its study, including a numerical simulations. 
First, two-clouds isothermal collisions are studied: for identical rotating clouds non-rotating 



clouds with different masses [141 ] and collisions of identical clouds with taken into account cool- 
ing/heating processes 15. [lit. 

Simulations with using great particle numbers (~10 5 — 10 6 ) and, therefore, with hight resolution 
require sufficient CPU consumptions even for supercomputers. In consequence of that, they carr yin g 
out either without integrating of energ y eq uation, but with barotropic state equation p = p(p) (l 71 ] . 
or simply in isothermal approximation [181 ] . 



In Monaghan & Varnas 15J] and Bhattal et al. [17j] the system of 48 and 1000 clouds corre- 
spondingly was considered. But even in such rich complexes the main role play just a pair collisions 

The most of such works were carried out without quantitative analysis of fragmentation. Thus, 



only in Bhattal et al. [13] and Gittins et al. 18[] the masses of fragments are evaluated. Additionally, 
in the last work qualitative analysis of IMF is given for which the fragments are defined as a single 
sink particles without any internal structure, etc [see QjJ, for more details]. 

In our work we carried out simulation of isothermal supersonic collisions of two identical molecu- 
lar clumps with aim of more detailed fragmentation study and mass distribution analysis of obtained 
fragments. For fragments finding we apply the modified variant of well-known " friends-of- friends" 
algorithm (FOF). 
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2 The Method 



2.1 Hydrodynamics 



For gas process modeling we use the Smoothed Particle Hydrodynamics method (SPH), which was 
independently suggested by Lucy 20] and Gingold & Monaghan (2l| . 

As in the well-known N-body algorithm, modeling system in SPH is represented by a set of 
computational particles which characterized by mass, position, linear and angular momentums and 
energy. Besides, each particle is diffusive one and additionally have some hydro- and thermodynam- 
ical parameters (such as density and thermal energy), defined for its "center" by interpolating from 
nearest neighborhood. Partially overlapping, SPH-particles gives the average weighing contributes 
to local properties of medium, according with their internal distribution profile. 

So, the smoothing estimate of any physical parameter /(r) in SPH is: 
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where the integral is taken over entire space and summation is over all N particles, h— so-called 
smoothing length, W(r — r', h) - interpolating smoothing kernel, which must satisfy to two simple 
conditions: J W(r — r', h) dr' = 1 and lim^^o W(r — r', h) = 5(r — r'). 

Yet another SPH features is that it need not in finite-differences for gradients calculating. Instead 
of this we must only differentiate the smoothing kernel and then smooth interesting parameter as 
ID: 
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The concrete form of kernel is chosen from reason, that non-zero contribution to local hydro- 
dynamical properties must take only the nearest particles. Smoothing length h define a scale of 
smoothing and so the spatial resolution. 

The great flexibility of SPH is provided by absence of any restriction neither on system geometry, 
nor on admissible range of varying of dynamics parameters, easiness of kernel change (that is equal 
to change of finite-difference scheme in grid methods), variety of equation symmetrization means, 
artificial viscosity forms, etc. 

Being hydrodynamics extension of earlier N-body method, SPH can naturally include it, as a one 
of its parts, for modeling collisionless components (as a dark matter). For handle with selfgravity 
field SPH use either grid methods to solve of Poisson's equation, or TREE-algorithm or, at 
last, direct summation of all particle-particle interactions by special hardware GRAPE. Also, SPH 
may be easy combined with other optional algorithms such as starformation, feedback etc [see for 
example El, 0, 0, • 

For more details see the goo d reviews by Monaghan 28], Berczik &z Kolesnik [2^ ]. Hernquist fc 



Katz [22f 
Berczik 
Miiller 



Hiotelis & Voglis [30J, and other works such as Lombardi et al. [31(, Thaker et al. 32|, 
2$ . Berczik & Kolesnik 33J], Carraro et al. [23], Monaghan & Lattanzio [34[, Steinmetz & 



J and of course Springel et al. [36( . 
We used the parallel version of freely available SPH code GADGET version 1.10 (GAlaxies with 
Y)ark matter and Gas infEracT) [i^], that destinate for numerical SPH/N-body simulating both of 
isolated self-gravitating system and cosmological modeling in comoving coordinates (with/without 
periodical boundary conditions). The GADGET v 1.1 distributive contain two codes: serial one, for 
workstations, and parallel - for supercomputers with distributed memory or PC-clusters with shared 



1 http : //www .mpa-garching.mpg. de/gadget 
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memory. Code is written on ANSI C and use standard MPI library (Message Passing Interface^ 
for parallelization. 

GADGET allow up to four types of N-body particles for collisionless fluids and one type of SPH 
particles for gas in modeled system. The main code features are spline-based smoothing kernel 



341 ]. which is second order accuracy (in since that </(r)>= /(r) + 0(h )), shear-reduced form of 



generally accepted Monaghan artificial viscosity 3l|, |32|, |28[ , Burnes & Hut (BH) TREE-algorithm 



for gravity calculation, leapfrog scheme as a time integrator [3CJ, l29|] and so. 

The code is completely adaptive in time as well as in space, owing to algorithms of individual 
time-steps and smoothing length for each particle. The last feature signify, that the number of 
nearest smoothing neighbors, which must influence on local medium properties is just constant in 
serial code version and nearly constant (in small enough bands in a few percents) in parallel version 

The time-step for each active particle calculated as 

Atx = (3) 
l a «l 

Additionally, there is yet another hydrodynamical criterion specially for SPH particles: 

K |(Vv)j| + max(cj, |vj|) • (1 + 0.6a visc ) ' 
and taken the minimal from both, (here a is acceleration, atoii ~ dimensional accuracy factor [see 



361 ] . etc ~ analog of Courant factor, a v i sc - artificial viscosity factor, c - sound speed). Besides, 
time-steps are bounded by their lower and upper limits: Ai m j„ and At max . 

For gravity calculation in GADGET were employed the standard BH oct-TREE, which is built 
for each particles type separately. For SPH particles it used also for neighbors search. Gravity 
potential expanded up to quadrupole order. As a cell-opening criterion may be used either standard 



BH-criterion 37|, or original one, based on approximative estimate of force error, caused by multiple 



expansion oneself [3 
2.2 Fragments finding 

As our purpose was to study the fragmentation in colliding clumps, we has not only to simulate 
the collision, but also to find an appropriate method for the determination of the extent of medium 
fragmentation and for the statistical analysis of fragments properties. 

We decided in favor of the well-accepted " friends-of- friends" algorithm (hereafter FOF). It was 
initially proposed for the separation of galactic groups and clusters in the spectroscopic redshifts 
surveys |38| . This simple and versatile algorithm is still widely used [3^ . l4o| . As original FOF 
method handled with points on coelosphere and distances between them it can be easily adapted 



for analyzing the results of numerical N-body/SPH experiments [4ll. l42f]. 

For each particle, that has not yet been subjected to the analysis procedure, all its neighbors 
located inside the sphere of a certain so-called "search-linkage" radius 5 are found. When there are 
no such neighbors, the particle is write off as isolated one (or as a background particle), otherwise 
this procedure is recursively repeat to all its neighbors are found, neighbors of neighbors and so 
on, until the whole list is exhausted. Then all the particles thus found are included in a cluster 
(fragment) if their amount is not below some threshold N m i n . In other words, if any two particles 
are separated by less then 5, they both belong to the same group (but only in the case when the whole 
group will collected no less then iV m j n members, else all these particles are reckoned as belonging to 
the background). 



2 http : //www-unix .mcs . anl .gov/mpi/ 
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The advantages of the FOF algorithm are its coordinate-free nature and so absence of any meshes 
and independence from problem geometry; no a priori assumptions are made as to properties of 
the sought-for groups of particles (such as shape, symmetry, density profile etc.), which is very 
convenient in the work with 3D models. 

The algorithm has only two "adjustment parameters": 5 and N m i n . The former parameter 
control the compactness and extension of the sought-for fragments, while the later one rejects random 
small groups of closely located particles like pairs and triples (the statistics suggests, that they 
comprise about 80 — 85 % of all groups found). Generally spiking, N m i n should be as small as 
possible, but not smaller then the number Nb of smoothing neighbors. 

Choosing of "search-linkage radius" (SLR) 5 is not all the clear, but obviously it should not be 
too large. It often taken in terms of the mean interparticles distances (if that's case the geometric 
mean is better, than arithmetic mean, because the existence of very distant, very compact groups 
can result in unjustifiably large 5). In this case, however, at different moments the fragments are 
separated with the use of different scales, i.e., 5 will depend from time. In addition, any mutual 
fragments moving off will influence the results of analysis. So, it is necessary to use the most stable 
criterion (as regards to interparticles distances varying) for SLR choosing, or set it by hand. 

Apart from all, the simple FOF method have some drawbacks. For example two groups which 
linked by a small thread will be identified as a single structure; or when there are several groups 
with strong different compactness and/or with internal fragmentation, the algorithm may either do 
not recognize the internal structure of some fragments, or it can let go some less compact groups at 
smaller 5. 

We modified the standard FOF algorithm to avoid such situations, by manual setting up of a few 
fixed SLR 5i and multistep fragments searching. All <5j are consequently applied by they increasing 
to nonidentified particles (i.e. writing off from previous step). So, first the most compact fragments 
will be detected with the smallest <5; next the less compact fragments are found with the new SLR 
among remained particles and so on. 

We took four radii 5 = 0.001,0.005,0.01,0.05 pc. The largest radius is approximately equal 
to the gravitational smoothing length (see below), and smallest radius was chosen for the densest 
fragments recognizing as a separate ones. The threshold particles number N m i n was equal to number 
of smoothing neighbors Nb = 40. 

3 Scaling and initial model 

GADGET allow to use arbitrary user's system of units, including dimensionless variables. As a base 
units, are units of mass, length and velocity. We choose 

[m] = lSEJIo, [/] = 1 pc, [v] = 1 km/s. 

Following this, unit of time is defined as [t] = [^[v]' 1 « 0.978 Myr, density - [p] = [M][l}~ 3 sa 
6.77-10 -20 kg/cm 3 (that correspond to number density 17.6 cm -3 , at choosen mean molecular weight 
- see below). The temperature measured in Kelvins, but everywhere in code use the internal specific 
energy u instead of temperature, which is measured in [v] 2 . In chosen units gravitation constant is 
G = 4.3 -10" 3 VJIq 1 • pc • (km/s), gas constant - 9t = 8.314 -10" 3 9K - (km/s) 2 - mol" 1 - K" 1 . 

We consider four models of supersonic collisions of two identical molecular clumps with various 
impact parameters j3 = b/R = 0, 0.2, 0.5, 0.75 (were b is linear shift of clumps centers, R - their 
radius) and therefore with different initial angular moment. 

Firstly, clumps' centers positioned at {— R; §;0} and {R; — 1;0}. There are no any internal gas 
motions in both clumps initially, but they moves to each other along A-axis each with velocity v 
(see Fig. [T]for example). 



5 



The mass 9Jt in both clumps is distributed inside cut-off radius R according to density profile 

We choose the follow basic parameters for each clump: 

9K = 2OOO9Jt , i? = 3pc, T = 20 K, u = 5km/s, fx = 2.3 
That is in good agreement with observational data about ISM in our Galaxy 0, 43|, 44, 0|. 



Additionally, each clump have further derived characteristics: 

► free-fall time: r// « 1.92 Myr, 

► crossing time: T cross = « 0.59 Myr, 

► isothermal sound speed: c « 0.27 km/s, 

► mean number density (at given /i = 2.3): n ~ 300 cm -3 , 

► ratio of thermal energy to gravitational one: Ex/\Eq\ ~ 0.09, 

► Jeans' mass: 971 j « 1 9JTq in center and « 10 971© at edge. 

For generation of initial configuration we divide up each clump into 100 shells of equal mass 
with boundary radii obeying to density profile ([5]). Next, we homogeneous and randomly distribute 
all particles in total number N among this shells in corresponding shares independently for each 
clump. It is need to note, that then any ideal symmetry was wittingly excluded and there was slight 
departure in positions of clumps centers from theoretical ones. 

The two set of calculations were carried out: one with total number of particles iV = 2x4000 and 
second - with N = 2x 8000. In all calculations we used the fixed number of smoothing neighbors 
N B = 40, a visc = 0.75, a toll = 1.00 km/s, a c = 0.01, At min = 10" 6 [t], At max = 0.75 [t] and 
gravitation softening e = 0.03 [see 13(1] . 

For calculation of all models we had used the supercomputing facilities of Hungarian National 
Supercomputing Centeid (SUN E10k/E15k supercomputers). 

4 The Test Runs 

Besides main simulations, we first fulfilled a few simple test-runs for study of code characteristics 
oneself. 

First test was adiabatic collapse of cold cloud with density profile ([5]) and homogeneous specific 
internal energy distribution 

K 

in natural dimensionless units system G = 9Jt= R = l. 

This test, named as " Evrard-test" , was sug gested in Evrard [IB] and now became a standard 



tool for SPH-codes verification [221 . |32j, l35l . [27j, |36|. It was calculated out with different particles 
number iV = (1, 2, 4, 8, 16, 32) • 10 3 till final moment nearly equal to 3 Tff. 

In the full agreement with cited works, at t ~ 0.7 — 1 Tff central bounce was occurred and during 
outwards shock wave going, significant dissipation of kinetic energy into heat taken place. Nearly at 
2.5 Tff in whole system the virial equilibrium between thermal and potential energies Et ~ —Eq)2 
had been set. Also, test showed the excellent energy conservation \8E/E\ < 0.1% even at small N. 

Generally speaking, the total CPU-consumption must determine by selfgravity calculation. Our 
results showed the typical for TREE-algorithms asymptotic of time consumption ~ 0(N log(N)), 
but the hydro dynamical consumption is even greater then for gravity. 



J http : //www . iif .hu/ szuper 
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As the next test, we calculated out the spherical isothermal collapse of analogous cloud with 
mass dJt = 2000 9Hq, radius R = 3 pc, but with temperature T = 20 K, as in main simulations. 
Particles number was set to N = 4000. After 1.56 Myr, that correspond to 0.81 Tff, central dense 
core with size of ~0.03 pc was formed close to centeiQ with thin extensive shell around, that have 



self-similar density distribution p oc r [see |46[] . 



The third test was fulfilled for check-up of comparative role of each time-step criterion ([3]) and 
and for choice of available accuracy factors. We calculate adiabatic collapse from " Evrard-test" 
with 32000 particles once with fixed a toU = 0m(GWl/ R) 1 / 2 and various a c = 0.1, 0.05, 0.01, 0.001, 
and another time - the wrong way, with fixed ac = 0.01 and various ettoii = 0.02, 0.05, 0.07. With 
changing of Courant factor ac the CPU-consumptions was increasing pro rata to its decreasing, 
that point to determinant role of criterion (QJ). Since energy error with ac decreasing approach 
to some limit (~ 0.05%), we choose - as a available from ratio accuracy/consumption - its value 
ac = 0.01. At such choice, variation of atoll even does not any affect neither on accuracy, nor on 
time-consumption. 

And in the end, we calculate one of clumps collisions - with iV = 2x4000, v = 5 km/s, (3 = 0.5 
during 1.3 Myr on 12, 4, 8, 16 and 32 SUN UltraSPARC III (1050 MHz) processors for verify of 
code parallelization.Lj (see tab. [[]) 

In general case, any parallel code may be split to so-called consecutive and parallel parts. There- 
fore it is easy to characterize a quality of code parallelization by relative share of this part in total 
time-consumption. For example, on parallel part fall 

1 - 1/k 



1 - 1/JVp 



where N p - processors number, fe > 1 — ratio of durations of task work on single CPU and on N p 
processors. For GADGET this value is greater then 93%. Also, we can evaluate the effectiveness of 
CPU loading, as e = k/N p (tab. 



5 The Results 

Head-on collision (f3 = 0) Because of extensively CPU-consumption this model was traced 
till 2.885 Myr only. However, already at 1 Myr for N = 2x8000 and at 1.76 Myr for TV = 2x4000 
the number density of central object amounted of 4.5 -10 10 cm" 3 and exceeded the "opacity limit" 
(10~ 13 g cm -3 = 2.6-10 10 cm -3 , see Masunaga & Inutsuka [43]). So, there was no reason to continue 
simulation in isothermal approximation at last for this object. 

The column-density images in XY plane and particles number densities vs. X-coordinate on 
two moment are shown on Fig. [2] 

Since the collision was going on between the denser clumps centers, the density in shock front had 
been strong increased and rapid fragmentation there was occurred, that yield the first maximum of 
fragments amount (see Fig. [6j). Next, this fragments began just as rapid merging one with another, 
giving a few massive central objects (minimum of fragments amount on the same plot). Meanwhile, 
the rest gas accreting on them, had been compressed and began intensively fragment too, that give 
the second peak of fragments amount. 

Finally, most of gas (> 70%) fell onto the flatten dense compact core of mass 2921.5 SDTq and 
number density up to 6.7 • 10 11 cm -3 . There are two non-coplanar discs with masses 537 and 
431.5 9JIq around this core. The external ring have three small subfragments (< 20 9#o) and entire 
system are embedded in thin extensive shell with mass 20 9JIq. 

4 due to character of initial configuration building, there is no exact coincidence with geometrical center (see sec. 
for more details.) 

5 for single CPU calculation the serial GADGET version was used. 
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The fact of forming of rotational system in head-on case may be explained by features of initial 
model building (see sec. [3j) and therefore there is a small impact parameter close to 0.01 pc. 

In calculations with half particles number the main picture remain similar to above, but the 
central body had a one single semi-destroyed most massive accretion disk abounded it. This "disk" 
also have one subfragment and entire system was embedded in thin shell. 

Collision with (3 = 0.2 At such small impact parameter clumps' centers, even if not col- 
lided, but passed close enough one from another and collision finished, as above, by forming of one 
massive zone of fragmentation^ in center (the regions in which fragments are formed by groups) and 
extent, mainly unfragmented diffusive big two-arm spiral with spread in ~30 pc and thickness from 
~ 3 pc in center to ~ 10 pc on periphery. In contrast to center, each arm consist from gas of its 
"own" clump. 

The some moment of evolution are illustrated in Fig. [3l But, at the end of simulated we have 
the massive flatten core with mass in 1141 and number density < 1.7 -10 10 cm" 3 in center and 
its disc-shell by mass 1763.5 9Jt©, which oneself have three petty subfragments (10 — 15 9Jt ). They 
are surrounded by internal two-armed spiral with size close to 1.5 — 2 pc, that have four fragments 
(164.25, 24, 16.75 and 12 9JIq). Rest fragments with masses nearly to 20 9JIq are detected in both 
arms of great spiral. 

In simulations with 8000 particles the external spiral was practically unfragmented with the 
exception of one fragment with mass 22.5 9JTq. 

The core looks like a one-armed spiral of ~ 0.75 pc and three satellites with masses in a tens 
solar masses. 

Collision with (3 = 0.5 At (3 = 0.50 central clumps' regions, as before, passed through 
shock front, has been formed and where at ~0.6 Myr appeared the first fragments with good mixed 
gas. Rest material, mainly from outer layers of clumps, go past far from "place of collision" and form 
large two-armed spiral structure about of 40 pc. The next fragments generation had been formed 
in outer sided of both arm after ~ 1 Myr since extensively fragments formation in shock front. 

Now collision gave 4-5 big zones of fragmentation, instead of only one such zone in two previous 
cases. They are the close pair in center, one single fragment about them and two less massive 
condensations on the both ends of spiral arms. 

The column-density images and X vs. Lg(n) plots at 0.83 and 1.47 ; Myr are shown on Fig. [H 

The central pair is two massive flatten cores of masses 563.25 and 475.5 OJT© with number 
densities beyond 10 9 cm" 3 , which are surrounded by accretion discs of masses 563.75 and 610 9Jt©, 
correspondingly and by order less dense. But, the first core yet have low massive thin flat shell 
(23.25 9JTq) around, and second one continuously (without space) transform to their disc-shell. 

In vicinity of central pair there is a third smaller dense fragment with mass 351.5 971© and 
number density < 2.7 -10 9 cm -3 . 

There are core 399.75 9Jt©, (n < 1.5 -10 9 cm -3 ) with double-shell of 73.75 9Jt© on the end of 
one arm, and another one 382.75 9Jt©, (n < 8.9 -10 8 cm -3 ) with shell in 54.5 SQt© on the other side. 

As for calculations with half mass resolution, in the first place spiral arms are mainly non- 
fragmented, also, both central core have accretion discs and surrounded by thin flat shells. 

Collision with (3 = 0.75 (see fig. illustration) This is largest impact distance we was 
considered. In difference from all previous models, now clumps' centers passed edgeways through 
forming shock front. Therefore, at chosen clumps velocities and such wide impact parameter the 
whole collision yield the greatest fragments amount, but with lesser degree of gas mixing since main 
shaking take place just in shock front. 

As before, first intense formation of fragments appeared at ~0.65 Myr in the shock front and the 
second fragments generation had been forming in outer side of both lengthy arms ~0.9 Myr later. 
But soon, almost all fragments and other gas go far from place of their "birth", so that the nearest 

6 "zone of fragmentation" or "condensing zone" 
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to center fragments, which are found on about 6 pc from center, are relatively massive dense disc 
with 9JT = 414 SDt©, n < 1.1 -10 9 cm -3 and smaller core in 51.5 DJIq with shell. The rest fragments 
have larger remoteness of 16 — 17 pc from center. 

On both ends of extent structure are concentrated by about 1500 9710 in fragments. There are 
core in 759.25 SDt© {n < 4.8 • 10 9 cm" 3 ) with small shell of 54.15 SOT©. On the another end we 
have pair of massive fragments groups too: two dense cores with masses 188.25 and 614 9Jl & with 
shells in 591 and 71.75 9JIq correspondingly. Besides there are yet a few smaller fragments here 
with masses less then 40 SO?© scattered around. 

In simulations with worst mass resolution we have more poor fragmentation, but in other respects 
after all, there are lesser qualitative distinctions between simulations results for N = 2 x 4000 and 
N = 2 x 8000, then in three previous models. 



6 Discussion. 

At chosen clumps parameters, their masses significantly exceed Jeans' mass, hence they are grav- 
itationally unstable. Theoretically each clump, taken separately, will collapse out during free-fall 
time-scale and all fragments, formed under contraction, in the end will have fallen onto the central 
object. But picture had been radically changed, if both clumps collide rapid enough one with an- 
other, so that collision time is lesser then free-fall time (in our models t ctoss : tj-j- = 1:3). Now, the 
evolution of forming fragments will depend from initial impact parameter /3 or (that is same) from 
initial angular moment of whole system. 

In spite of (3 varying in a wide range (and independently from N) in all models began intensive 
gas fragmentation after nearly t^-oss from start of collision, that is in good agreement with such 



previous work as Gittins et al. [18[ and Bhattal et al. [17]]. The number density at this time never 



exceed opaque limit in ~2.6 -10 cm (see sec. H]); hence the isothermal approximation is quite 



available, at last, during this period [see l46l. |47| 



Some general information about gas fragmentation is given in tab. [5J But it is necessary to 
note here, that owing to improvement of resolution, with invariable SLR and N m i n (see sec. 12.21) 
at greater particles number N fragments will be detect earlier and in a larger amounts, then in 
models with the worst resolution. Therefore, in the tab. [2] are given moments of begin intensive 
fragmentation in frontier shock front forming during collision, instead of time of first fragments 
detection. 

Since we have carried out comparatively low-resolution modeling we concentrate in this article 
mainly on integral characteristics of system. The main attention we devoted to the general properties 
of fragmentation - mass spectra of condensation, their mass share in whole system, etc. 

There is a clear correlation between impact parameter (3 and quote of remained fragments to its 
maximal quantity during evolution: k = N final /^max oc (3 (see tab. [2]). This hint at interclumps 
collisions may play role not only as a trigger mechanism in ISM fragmentation/starformation, but 
also as stabilizer of formed fragments system and set the rate of this processes. 

We was not built simple, traditional differential mass function, since they showed most sensi- 
tiveness to chosen mass interval, quantity of fragments and so, then cumulative one. 

The integral mass spectrum can be derived from differential one by its integration either in 
bounds (0, m), that give the cumulative mass function iV(< m), i.e. the quote of fragments with 
masses less m, or in bounds (m, oo), give the inverse cumulative function iV(> m) - the quote of 

dN 
dm 

dm 



fragments with masses greater then m. If differential distribution ^ is approached by generally 
accepted power law ^ = km~ a (k,a= Const, a>0), then we have 



m 

N (< m ) = I ~77~d£ = T~~ — ml ~ a + c, 
J d£ 1 - a 

o 
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(where c = is constant, connected with minimal fragments mass) and 



oo 

dN _ k x _ a 



N(> m)= I -—d£ = -- m 

i d( 1-a 

m 

where we had taken into account, that iV(> m) is a decreasing functiorjj|. 

Since we have not a many data points in our spectra, we choose the inverse cumulative function 
iV(> m) = 1 — iV(< m) oc m~ 7 (7 > 0) for data fitting, as it have only two free parameters, instead 
three ones in another variant and may be linearized by transform to Zog-variables before fitting. The 
obtained mass function are given on Fig. [7J and their slopes are summarized in tab. [20 

The trace of fragments amount and total fragmentation index rj (that we defined as ratio of all 
fragments mass to mass of whole system: rj = 9JT/-/9JT) for all evolutionary sequences are given on 
Fig. [6] together with slopes 7 of fragments mass spectra. Above all, it is clearly seen the great 
similarity of curves at the same (3 for both N, with the exception of differences in fragments amount 
on behalf of better resolution and limit values of 7. Also, it is seen, that always fragments amounts 
firstly strongly increase during about 0.5 — 0.6 Myr at J3 = 0, 0.2 and ~ 1.5 Myr at (3 = 0.5, 0.75 up 
to its maximum but next begin decrease, moreover at small (3 = 0.00, 0.20 most rapid, then in two 
other cases. As to the total fragmentation 77, it rapid increase together with fragments amount, and 
next its grow practically stopped, not changing to decreasing. 

This fact may be explained as that intensive gas disintegration begin at nearly T cross and next 
had been broken out after a time about 1.5 — 2 Myr. But existing fragments continue stick one 
with another and accrete surrounding gas. Apart from this, the great fragmentation was being 
achieved at larger impact parameters, with one exception of head-on collision, that gave mainly 
total fragmentation. 

Finally, we had built the dependence of cumulative mass functions slope us. time for each calcu- 
lated model, that is present on Fig. \U\ by solid line. 

As it seen, 7 smooth rush to certain limit, that is —0.5... — 0.6 in models with N = 2 x 4000 and 
close to -0.4 in models with N = 2 x 8000. 

But, since 7 is continual vary and admissible of our isothermal approximation is restricted in 
time, we must dwell on some "typical" its values for represents of our results and comparisons of 
mass spectra slopes between different models. It is naturally to tie up the typical 7 with certain 
typical time moments in connection with process of fragmentation. We chose two such benchmark 
moments - T0.8, when 80%-fragmentation is reached, and 2ro.s (Fig. UJ). 

When 80% of gas got to fragments the mass spectra slopes are nearly 0.5 for first two models 
with (3 = 0, 0.2 and 0.7 for two other models (with estimate error less then ±0.03). After yet 
To. 8 period the difference between slopes for collisions with various (3 had been rub off, all spectra 
became flatten and amount to ~ 0.4 — 0.6. 

The resulting slopes of mass spectra are in good agreement with many observational data for 
molecular clouds and clumps in our Galaxy, especially for large collisions with (3 = 0.5, and 0.75. 
Thus, Casoli et al. [s| give 7 = 0.59 for 250 clumps with masses from 4-10 2 to 2 -10 5 in Perseus arm. 
Simon et al. [48| from FCRAO and Milky Way Galactic Ring Survey in 13 CO spectral line received 



7 = 0.80 ±0.10 for central Galaxy region; Heyer, Carpenter & Snell 49J] from FCRAO Outer Galaxy 
Survey in lines 12 CO/ 13 CO received 7 = 0.80 ± 0.03; Tatsematsu et al. |50|] from CS(l-O) emission 
for Orion A GMC gave slope 0.6 ± 0.3; Stutzki & Glisten [jj by analysis of high-resolution maps 
of M17 SW in C 18 0(2-l), C 34 S(2-1) and C 34 S(2-1) lines calculated 7 = 0.72 ± 0.15; Kramer et al. 



7 Of course really the integration going in finite mass interval from to m ma x 

8 However, it is necessary to note, that as we used multiscaled fragments searching, the massive fragments with 
complex internal structure had been split in a few lesser subfragments. Owing to this fact, the whole mass spectrum 
became flatten. 
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Ill ] for seven molecular clouds (also including M17 SW with the same slope) found mainly similar 
slopes 7 = 0.6 — 0.8 in a wide mass range 10~ 4 — 10 4 SOT©; 

In spite of not great number of models (in all four) they have divided onto two types with 
different characteristics - collisions with small and large impact parameter (0 — 0.2 and 0.5 — 0.75, 
correspondingly). In first group, owing to insufficient initial angular moment, there are single 
dominated by mass center of fragmentation with more or less developed satellites system, spirals etc 
have been formed. The fragments amount, rapid decrease, after its maximum, but mass distribution 
slopes are ~ 0.5 after ro.s and ~ 0.4 after twice period. 

Another group - collisions with large /3 - lead to form essentially more advanced set of fragments. 
Now there are a few nearly equivalent centers of fragmentation, so each of they have its own mass 
dominant component and several satellites. The fragments amount curves are more flatten and mass 
distribution are steepen: 0.7 — 0.8 and 0.5 — 0.6 on moments ro.s and 2tq.8 correspondingly. 



7 Conclusions. 

The main our results are follow. 

1. After nearly T cross after start of collision in all models begin intensive gas fragmentation. The 
fragment quantity rapid increase up to its maximal value during about 0.5 Myr at (3 = 0, 0.2 
and ~ 1.5 Myr at /3 = 0.5, 0.75. Next, existing fragments continue stick one with another, 
but mainly no more fragments formed. Owing to this fact fragments amount decrease (more 
smooth at lager (3), but their total mass share stay almost constant, from now. 

2. With the exception of head-on collision, the total gas fragmentation r\ increase proportional to 
initial impact parameter. But, at that time the maximal degree of fragmentation (> 98.6%) 
had been reached in head-on case. 

3. The quote of remained fragments to its maximal quantity during evolution k is direct pro- 
portional to impact parameter, what mean that clouds/clumps collisions may play role not 
only as a trigger mechanism in ISM fragmentation, but also as stabilizer of formed fragments 
system. 

4. Finally, obtained mass spectra of fragments, that have typical slopes 0.5 — 0.7 (see above) 
accord with many observational data about clouds and clumps in the Galaxy. The mass spectra 
became flatten with time, but steepens with impact parameter (in time range [ro.s; 2ro.s])- The 
difference between two type of collision had been rub off during second ro.s period (from ro.s 
to 2ro.8) and was set to 0.4 — 0.6. 
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N s tep 


ttotal 


tgrav 


tfiydro 


speed 


Tp 


e 


1 


177,482 


17,461.27 


6,455.20 


10,285.22 


10.164 






2 


167,882 


9,254.15 


3,778.87 


5,033.96 


18.141 


0.940 


0.943 


4 


164,704 


5,059.01 


1,903.46 


2,723.69 


32.557 


0.947 


0.863 


8 


168,956 


2,741.42 


925.20 


1,549.85 


61.631 


0.963 


0.796 


16 


167,245 


2,170.76 


630.79 


1,269.27 


77.044 


0.934 


0.503 


32 


169,125 


2,346.10 


570.27 


1,328.61 


72.088 


0.894 


0.233 



Table 1: CPU consumptions for test. N proc - CPUs number (SUN UltraSPARC III 1050 MHz), N st£p - time- 
steps number, ttotal— the total time consumption on whole task in seconds, t grav - time consumptions on gravity 
calculation (including tree-building) in seconds, t hydro - the same on hydrodynamics (including neighbors search), 
speed - the total calculation speed in sps (steps per second), r p - the share of parallel part of code in ttotal, e - the 
effectiveness of CPUs loading (see text). 
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/? 


begin fragmentation 


end of simulation 


7 


r | n max N f 931/ 


t n max \ Nf \ Mf T] | ft 


Tq.8 1 2r . 8 


N = 2x4000 


0.00 


0.59 


4.68 -I0 b 


1 


114.5 


2.89 


4.68 -10 12 


4 


3945.0 


0.986 


0.257 


0.55 


0.19 


0.20 


0.59 


5.20 -10 5 


2 


107.0 


5.87 


2.13 -10 1 " 


8 


3304.0 


0.826 


0.351 


0.37 


0.40 


0.50 


0.59 


4.94 -10 5 


1 


21.0 


5.87 


1.14. 10 11 


17 


3519.0 


0.880 


0.390 


0.79 


0.58 


0.75 


0.64 


1.04 -10 5 


1 


26.5 


5.87 


4.16-10" 


17 


3578.5 


0.895 


0.489 


0.84 


0.61 


N = 2x8000 


0.00 


0.54 


3.70 -10 7 


11 


244.50 


2.89 


6.68 10 u 


9 


3943.75 


0.986 


0.235 


0.55 


0.38 


0.20 


0.49 


1.68-10' 


3 


41.25 


5.87 


1.65 -It) 1 " 


13 


3451.50 


0.863 


0.381 


0.52 


0.43 


0.50 


0.54 


3.30-10' 


4 


59.75 


5.87 


8.98 -10 1 " 


16 


3600.50 


0.900 


0.548 


0.71 


0.39 


0.75 


0.59 


3.46 -10 7 


3 


48.50 


5.87 


9.44-10" 


23 


3698.25 


0.925 


0.486 


0.73 


0.50 



Table 2: The parameters of fragmentation for models with both N. r- time in Myr, n max - the maximal 
number density (in cm~ 3 ) at this moment, Nf- fragments amount in system, 371/ — mass of all fragments in SHq, 
r\- system fragmentation, k— the quote of remained fragments to its maximal quantity, to. 8- the moment, when 
80%-fragmentation is reached, 7- the integral mass function slope (dN/dm <x m -7 ). 
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Figure 2: Column-density images along Z-axis in large and small scales as well as common logarithms 
of number density vs. particles X-coordinate in model with (3 = 0. The top row is for time 
t = 0.78 Myr and the bottom one - for t = 1.03 Myr. 
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Figure 3: Column-density images with X — lg(ra) plots for /? = 0.2. The top row is 
t = 0.73 Myr and the bottom one - for t = 1.03 Myr. See notes on Fig. [2] 
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Figure 4: Column-density images with X — lg(n) plots for (3 = 0.5. The top row is for time 
t = 0.83 Myr and the bottom one - for t = 1.47 Myr. See notes on Fig. [2] 
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Figure 6: The fragments amount Nj (boxes), total fragmentation r/ (dotted line) and slope of frag- 
ments mass spectra 7 (solid line) vs. time for N = 2x 4000 ( top ) and N = 2 x 8000 (bottom ). The thin 
solid vertical line and two dash-dotted lines signs the Tff, tq.s and 2tq.8 moments correspondingly. 



22 



P= 0.00; 



t= 1.42 anil 2.84 Myr 



3= 0.00; 



t= 1.37 and 2.74 Myr 




Figure 7: The inverse cumulative mass spectrum for models with N = 2 x 4000 (left) and N = 2 x 8000 
(right). Solid lines for time moment To. 8 and dotted - for 2tq.8 
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